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Abstract: We present forecasts for the accuracy of determining the parameters of a 
minimal cosmological model and the total neutrino mass based on combined mock data for 
a future Euclid-like galaxy survey and Planck. We consider two different galaxy surveys: a 
spectroscopic redshift survey and a cosmic shear survey. We make use of the Monte Carlo 
Markov Chains (MCMC) technique and assume two sets of theoretical errors. The first 
error is meant to account for uncertainties in the modelling of the effect of neutrinos on 
the non-linear galaxy power spectrum and we assume this error to be fully correlated in 
Fourier space. The second error is meant to parametrize the overall residual uncertainties 
in modelling the non-linear galaxy power spectrum at small scales, and is conservatively 
assumed to be uncorrelated and to increase with the ratio of a given scale to the scale of 
non-linearity. It hence increases with wavenumber and decreases with redshift. With these 
two assumptions for the errors and assuming further conservatively that the uncorrelated 
error rises above 2% at A; = 0.4/i/Mpc and z = 0.5, we find that a future Euclid-like 
cosmic shear/galaxy survey achieves a 1-a error on Mj^ close to 32 meV/25 meV, sufficient 
for detecting the total neutrino mass with good significance. If the residual uncorrelated 
errors indeed rises rapidly towards smaller scales in the non-linear regime as we have 
assumed here then the data on non-linear scales does not increase the sensitivity to the 
total neutrino mass. Assuming instead a ten times smaller theoretical error with the same 
scale dependence, the error on the total neutrino mass decreases moderately from a[M^) 
= 18 meV to 14 meV when mildly non-linear scales with 0.1 /i/Mpc < k < 0.6/i/Mpc are 
included in the analysis of the galaxy survey data. 
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1. Motivations 

Several ambitious ground-based and space-based galaxy surveys have been planned for the 
next decade (e.g. SKA^, LSST^), or are about to take place (e.g. DES^). One of the most 
ambitious approved missions, the Euclid^ satellite [1], is expected to be launched by ESA in 
2019. It will combine a galaxy redshift survey with weak lensing observations, measuring 
the matter power spectrum and the growth of structure with unprecedented accuracy. 
This will offer a unique opportunity to improve measurements of cosmological parameters, 
including the neutrino mass, known to slow down structure formation on intermediate and 
small scales [2], as well as constraints on dark energy and modified gravity models. 

Recent constraints on the total neutrino mass appear to have converged on an upper 
limit of about 0.3 eV at the 95% confidence level (e.g. [3, 4, 5, 6, 7, 8, 9]), with the 

"^http:/ /www. skatelescope.org/ 
^http:/ /www. lsst.org/lsst/ 
^http: / /www. darkenergysurvey.org/ 
^http:/ /www. euclid-ec.org 
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notable exception of Lyman— a forest data, which gives an even lower bound of O.lTeV 
[10]. These constraints rely on a combination of data from Cosmic Microwave Background 
(CMB) experiments such as WMAP, Baryonic Acoustic Oscillations (BAOs), SuperNovae 
(SN) distance moduli, galaxy clustering and cosmic shear (especially from the SDSS^ 
and CFHTLS^ surveys). Data sets provided by Large Scale Structure (LSS) are partic- 
ularly important, since they are able to probe scales and redshifts affected by neutrino 
free streaming both in the linear and non-linear regimes. Neutrino oscillation experiments 
provide a lower bound of 0.05eV on the total neutrino mass, meaning that the allowed 
range is now significantly squeezed by cosmological data, and well within reach of future 
planned surveys. 

Several forecasts have already been published on the sensitivity of Euclid to cosmo- 
logical parameters, with a focus on dark energy, modified gravity, the neutrino mass, or 
other extensions of the minimal ACDM model (see e.g. [11, 12, 13, 14, 15, 16, 17, 18, 19]). 
However reliable forecasts are difficult to obtain; interpreting Euclid data on small (non- 
linear) scales will require a more accurate modeling of systematic effects than is currently 
achievable. This is true for both non-linear corrections to the matter power spectrum, and 
for effects specific to each survey. In the case of the galaxy redshift survey, for instance, 
redshift space distortions and scale-dependent bias. In the case of the cosmic shear sur- 
vey, noise bias in shape measurements [20]. Some authors have pointed out that without 
considerable progress in modeling these effects, the sensitivity to cosmological parameters 
might degrade considerably (see e.g. [18]). 

Current forecasts tend either to incorporate only linear scales and neglect these sys- 
tematics, or to include a small range of mildly non-linear scales and model systematics by 
including nuisance parameters which are then marginalized over. Introducing such nuisance 
parameters (for instance, in order to describe redshift-space distortions) still assumes that 
we can predict the shape of these effects, and reduce them to a simple family of curves. 
Hence, this approach is not the most conservative. 

On top of this, many forecasts are affected by a methodology issue: apart from two 
recent works [19, 21], they are based on a Fisher matrix technique, whose results depend 
on the step chosen in the calculation of numerical derivatives of the spectrum with respect 
to the parameters (see e.g. [21, 22]). 

The present forecast has three objectives: 

• First, we wish to use a reliable forecast method for the sensitivity of a Euclid-like 
survey to ACDM parameters and to the total neutrino mass, based not on Fisher 
matrices, but on a parameter extraction from mock data with Markov Chain Monte 
Carlo (MCMC). This goal has also been achieved very recently by [19], although with 
a different approach for modeling the galaxy redshift survey. To our knowledge, the 
present analysis is the first MCMC forecast of a Euclid-like galaxy redshift survey 
using as an observable the power spectrum P{k) in wavenumber space. 

^http:/ /www. sdss.org/ 

^http:/ /www. cfht.hawaii.edu/Science/CFHLS/ 
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• Second, we wish to incorporate non-linear corrections using the most accurate avail- 
able fitting formula accounting for neutrino mass effects, namely the version of 
HALOFIT [23] presented in Ref. [24]. This formula has been obtained by fitting to a 
suite of N-body simulations which incorporate neutrinos as free-streaming dark mat- 
ter particles, using the code first presented in Ref. [25]. The error in this formula 
specific to the neutrino mass was estimated by Ref. [24] to be Gaussian, with squared 
variance 

where fjy = uj^jura and kcj{z) is the non-linear wavenumber as defined and computed 
in HALOFIT. We include this in the likelihood as a fully correlated error, as described 
in detail in Appendix A, associated to a unique nuisance parameter. 

• In order to obtain conservative results while keeping the analysis simple, we will com- 
bine this correlated error with a second uncorrelated error. This second uncorrelated 
error is assumed to account for extra uncertainties in our approximate modeling of 
non-linear corrections, redshift space distortions, scale-dependent bias and other sys- 
tematic effects. By assuming an uncorrelated error on each data point, we remain 
more conservative than if we marginalized over a small set of nuisance parameters 
representing several types of fully correlated errors. Throughout this work, we as- 
sumed for convenience that the relative theoretical error on the power spectrum was 
given by Eq. (1.1), with fy replaced by a constant factor, by default 0.05. This error 
grows smoothly from zero on linear scales up to 5% on deeply non-linear scales. For 
a concordance cosmology and at redshift z = 0.5, it reaches 1% near A: = 0.1 /iMpc~^ 
and 2.3% around k — 0.6/iMpc~^. We assume that ten years from now, this will 
provide a reasonable description of the total uncertainty coming from all systematic 
effects in each of the two surveys. Occasionally, we will consider the effect of dividing 
the magnitude of the error by two or ten, to evaluate the effect of better control of 
non-linear systematics. We emphasise that the exact form of the uncorrelated error is 
obviously just an educated guess and that a different k-dependence will e.g. influence 
the assessment of how useful pushing to smaller scales will be. Of course, introducing 
a fully uncorrelated error (or alternatively, form filling functions as in [26]) is very 
conservative in that it assumes that no modeling of systematics is accurate enough. 
In several years from now, it might become realistic to model most systematics with 
several types of correlated errors, and to reduce the residual uncorrelated theoretical 
error to a smaller level than assumed in this work. 



2. Galaxy redshift survey 

Throughout this paper, our fiducial model is chosen to be a flat ACDM model with three 
degenerate massive neutrino species. The fiducial parameter values are taken to be cj^ = 
0.02258, ujc = 0.1109, As = 2.43 x 10"^ (pivot scale K = 0.05 /iMpc"^), Us = 0.963, 
h — 0.710, Zreio = 10.3, vrty — 0.07 eV (so M^, — 0.21 eV). For the power spectrum of 
the mock data, we could take directly the fiducial power spectrum, or generate a random 
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spectrum realization corresponding to the same model. As illustrated in [22], the two 
options lead to the same forecast errors, so for simplicity we assume an observed power 
spectrum equal to the theoretical power spectrum of the fiducial model. 




Figure 1: Observable spectrum (top) and relative error on this spectrum (bottom), for the first 
redshift bin (left) and last redshift bin (right) of a Euclid-like galaxy redshift survey The quan- 
tity displayed in the top is the galaxy power spectrum P^(/cref , /i, ^) as a function of the fiducial 
wavenumber /cref, for fixed redshift and perpendicularly to the line of sight (/i = 0), rescaled by the 
inverse squared bias b{z)~'^ and by a factor H (z) / D a{z)'^ : it is therefore a dimensionless quantity 
The upper plots show a comparison between a model with massless neutrinos and our fiducial model 
{Mjy = SrUiy = 0.21 eV). Solid lines are derived from the non-linear matter power spectrum using 
the updated halofit version of ref. [24] , while dashed lines are derived from the linear power spec- 
trum. The lower plots show the part of the relative error coming from observational or theoretical 
errors only (cosmic variance is included in the observational error). In these plots, the individual 
l-<j error on each data point has been rescaled by the square root of the number of points, in such 
a way that the edges of the error bands correspond to a shift between theory and observation lead- 
ing to Ax^ = 1, when only the observational or theoretical error is incorporated in the likelihood 
expression. In these lower plots, we also show for comparison the ratio between a massless model 
and a model with the minimum total mass allowed by neutrino experiments, Mj^ = 0.05 eV. 



We fit the mock and Euclid-like spectra using the MCMC code MontePython [27]. 
MontePython uses the Metropolis-Hastings algorithm like COSMOMC [28], but is in- 
terfaced with CLASS [29, 30] instead of CAMB [31], is written in python, and has extra 
functionality; it will soon be released publicly, including the Euclid-like likelihood codes 
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used in this work. 

Technical details of the assumed likelihood and our analysis are presented in Ap- 
pendix A. Let us summarize here the essential points. As in most of the recent Fisher- 
matrix-based forecasts, we assume that the reduced data is described by a set of observable 
power spectra P^^^{k^Qi^ fi^ z)^ related to the familiar non-linear matter power spectrum 
PNL{k, z) in a non-trivial way in order to take into account redshift space distortions, lin- 
ear light-to-mass bias, spectroscopic redshift errors and the Alcock-Paczynsky effect (see 
A.l). Of course, this modeling is imperfect: for this reason we introduce a theoretical 
error. For instance, we do not take into account galactic feedback [32], assuming that this 
contamination can be predicted by simulations up to the level of our residual theoretical 
error function. The arguments fcref and ji of the observable power spectrum stand respec- 
tively for the observed wavenumber assuming the fiducial cosmology, and the cosine of the 
angle between the observed wavevector and the line of sight. We assume sixteen redshift 
bins with mean redshift ranging from 0.5 to 2, and bin widths of Az = 0.1. For a fixed 
theoretical model, each observed value of P^^^ in a bin centered on the point {k^Qi^fi^z) 
follows, to a good approximation, a Gaussian distribution with variance 



^ ^ ^ref ^survey ^^Mi^^fcref/fcref J ^ ^ 

where dji is the size of the bins in fi space, and [dfcref/^ref] the size of the logarithmic bins 
in wavenumber space (see A. 2). The characteristics of the survey are encoded in V^urvey, 
the survey volume, and n^, the comoving number density of galaxies accounting for shot 
noise (see A. 3). Hence, if for every observed data point the theory and the observation 
differed by this amount, the effective would increase with respect to its minimum value 
by the number of data points, namely 

_ ^ 2 ln(fcniax/femin) 2) 

where B is the number of redshift bins. 

To illustrate this error, in figure 1, we show the relative error bar on the observed 
spectrum in the first and last redshift bin, assuming no additional theoretical error. For 
the purpose of comparing with the theoretical error introduced below, we do not show 
as usual the error corresponding to a one-sigma deviation for each given data point; we 
divided each error by \fN ^ in such a way that the edge of the error band corresponds to a 
deviation between the observed and theoretical spectrum leading to Ax^ = 1. Note that 
the displayed quantity ■^/S.P^^^ / [P^^^ \fN) does not depend on the width of the bins in 
(fcref, /i, z) space, but only on P*^, 1/^urvey and Ug. 

We incorporate the theoretical error in the likelihood in the way described in section 
A. 4. In few words, this error is normalized in such a way that a shift between theory 
and observations by a relative amount a (the quantity defined in eq. (1.1)) leads to an 
increase of the by one. This is achieved simply by adding a term N{aP^^)^ to the total 
error variance. Figure 1 shows the relative theoretical error on the observed spectrum, 
normalized in such a way that the edge of the error band corresponds to a deviation 
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between the observed and theoretical spectrum leading to Ax^ = 1 when the observational 
error is switched off. These edges are directly given by ±q:. 

We see in this figure that our assumption for a leads to an error of 1% at /c = 0.1/iMpc~^ 
and 2.5% at A: = 0.6/iMpc~^ for the first redshift bin centered on z — 0.5. For the last 
redshift bin in the galaxy survey, centered on z — 2, non-linear corrections appear on 
smaller scales, and the error is only 1% at A: = 0.6/iMpc~^. 
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Table 1: Marginalized l-a error for each model parameter, in a fit of Planck + Euclid-\^ke galaxy 
survey data. The different lines correspond to different choices of /cmax, to the inclusion or not of 
the global uncorrelated theoretical error (un. err.), divided by ten (l/io), by two (l/2), or full (•), 
to that of the specific neutrino-related correlated error (co. err.), and to the use of the non-linear 
or linear power spectrum. The models with correlated error have one more nuisance parameter 
not shown here, with unit 1-cr error. 

We performed several forecasts for a combination of Planck data and a Euclid-like 
galaxy redshift survey data. It should be stressed that the characteristics of Euclid are 
not yet finalized. Our choice for V^urvey ^.nd ng{z), detailed in A. 3, should be taken as 
indicative only. For Planck, we follow the method presented in [22] and do not include 
lensing extraction. For the experimental Planck sensitivity, we use the numbers presented 
in the Planck Bluebook^. This is a rather conservative model since the sensitivities are 
based on 14 months of observations instead of 30. 

The differences between our forecasts reside in the maximum wavenumber, equal to 
fcmax = 0.1 or 0.6/iMpc~^, and in various prescription for the theoretical error: no error 
at all, the uncorrelated error described above and in A. 4 (divided by ten, by two, or full), 
or additionally the correlated error accounting for neutrino-mass-related effects (described 
in A. 5). Since we are using an increasing theoretical error on non-linear scales, we expect 
the amount of information contained in the data to saturate above some value of kma^'- 
this is the reason we can consider such a high value as 0.6/iMpc~^. We did not try even 
higher values, first because our result would not change, and second because our forecast 

'^http://www.rssd.esa.int/SA/PLANCK/docs/Bluebook-ESA-SCI(2005)l_V2.pdf, page 4, Table 1.1 (us- 
ing only the best three HFI channels: 100, 143 and 217 GHz). 
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would become unrealistic: deep in the non-linear regime, the Gaussian assumption for the 
likelihood breaks down. 

Our results are presented in Table 1. Parameters like cj^ and Zreio are well determined 
by CMB data, and their forecast error depends very mildly on our different assumptions. 
For other parameters, the redshift survey plays a crucial role in removing parameter degen- 
eracies. In that case, even with fcmax = 0.1 /iMpc~^, including the uncorrelated theoretical 
error makes a difference: the parameter sensitivity degrades by up to 70% for h. The 68% 
neutrino mass error bar degrades by 40%, from cF{My) = 0.018 eV to a{My) = 0.025 eV. 

Assuming only this uncorrelated error, the cases fcmax^O.l /iMpc~^ and fcmax^O.G /iMpc" 
give almost the same results. Hence, our assumption for the theoretical error magnitude 
is such that most of the information is contained on linear scales. Thanks to realistic (or 
at least conservative) assumptions for the theoretical error, the results of our forecast are 
nearly independent of the cut-off fcmax- Without a theoretical error, increasing fcmax to 
0.6/iMpc~^ would lead to a spectacular (but totally unrealistic) decrease of the error bars, 
with cr(M^) = 0.0059 eV. 

If we are more optimistic and half the uncorrelated error, the error bars decrease 
marginally, as can be seen in the Table (lines starting with "1/2"). The error on the 
neutrino mass only decrease by ^ 10%. Assuming no error at all implies that the spectrum 
can be predicted up to the 0.1% level or better on small scales. In comparison, assuming 
a precision of one percent is not very different from assuming two percent. With the 
halved error, the sensitivity to the neutrino mass increases from a{M^) — 0.023eV to 
(j{My) — 0.022eV when including data in the range from 0.1 to 0.6/iMpc~^. 

Finally, in a very optimistic forecast with an error ten times smaller, we start to see 
how extra information can be extracted from non-linear scales; the error decreases from 
cr(M^) = 0.018 eV to cr(M^) = 0.014 eV when pushing fc^ax from 0.1 to O.G/iMpc"^. 

The inclusion of an additional correlated error accounting for neutrino-mass-related 
systematics has a negligible impact on our results. In our forecast, the uncorrelated and 
correlated part of the error have similar amplitudes and the same shape; however the 
uncorrelated error allows much more freedom and thus leads significantly more conservative 
results: this explains why the correlated error has a comparatively small effect. It should be 
stressed that our results depend not only on the assumed error amplitude at a given scale 
and redshift, but also on the wavenumber dependence of the error function a. Different 
assumptions, with a steeper or smoother step in the error function around the scale of non- 
linearity, would lead to different forecasts. In particular, as already mentioned the actual 
benefit from pushing to smaller, non-linear scales depends on the assumed k-dependence 
of the residual uncorrelated theoretical error. 

For the case with /cmax = 0.6 /iMpc~^ and no neutrino-related correlated error, we show 
the one and two-dimensional posterior probability on cosmological parameters in figure 2. 
We see several pronounced parameter degeneracies. For instance, the neutrino mass is very 
correlated with ujc and h. This suggests that further progress could be made by including 
extra data sets, such as direct measurements of the Hubble parameter, the cluster mass 
function, supernovae luminosity, 21-cm anisotropics, and so on. 

Our results are consistent with those of [14, 17], although a direct comparison is 
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Figure 2: Marginalized posteriors and two-dimensional probability contours in a fit of Planck 
plus a Euclid-like galaxy survey mock data, with k^^x = 0.6hMpc~^ and a global uncorrelated 
theoretical error (second line starting from the bottom in Table 1). 



difficult, since those authors include two extra parameters, wq and Wa, in their forecast. 
The results in Table 2.1 of [33], based on the same cosmological model, match our prediction 
in the case with no non-linear scales and no theoretical error included. A similar sensitivity 
was found by [19] for a Euclid-like photometric redshift survey, referred to as "eg" in their 
Table 2. However, this reference presents other results based on even more conservative 
assumptions than ours. We assumed that the bias function for each redshift bin could 
be determined in advance (up to corrections on non-linear scales contained in our global 
theoretical error). This assumption has also been made in most recent forecasts, since both 
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N-body simulations and higher-order statistics in the real data allow the prediction of the 
redshift-dependent bias of a given population of galaxies, at least on linear scales. Were 
this approach found to be unreliable, it would be necessary to marginalize over the linear 
bias in each redshift bin, b{zi). Ref. [19] did such a marginalization in the runs called "cgb" 
and "cgbl" , with no prior at all on each b{zi). They found roughly the same error bar on cjc 
and h than in our forecast with theoretical error, but a much larger error on the neutrino 
mass. However, it seems unlikely that at the time when Euclid data will be analyzed, no 
information at all will be available on the linear bias of the observed population of galaxies. 

3. Cosmic shear survey 

For the case of a Euclid-like cosmic shear survey, we stick to the same fiducial model and 
methodology as in the previous section. The likelihood is now a function of the observed 
lensing power spectrum Ci in harmonic space and for each pair ij of redshift bins, 
taking into account photometric redshift errors and shot noise (for details, see B.l and 
B.2). We assume experimental sensitivities summarized in B.3, and cut the observations in 
five redshift bins covering the range < z < 3.5 (although a negligible amount of galaxies 
contribute between 3 and 3.5). We do not take into account intrinsic alignment, assuming 
that this contamination can be removed up to the level of our residual theoretical error 
function [34, 35]. 

As explained in detail in B.4, there is a small technical difference between the likeli- 
hood of the galaxy survey and the shear survey in the way we incorporate the uncorrelated 
theoretical error. For the galaxy survey, the theoretical error was encoded as an extra con- 
tribution to the total error variance. This can be justified mathematically by marginalizing 
over one nuisance parameter for each data point. The shape of the galaxy survey likelihood 
allows for an analytical minimization over each nuisance parameter, in such a way that nui- 
sance parameters do not appear explicitly in the final likelihood. We found that no such 
scheme is accurate enough in the case of the (chi-square type) shear likelihood. Hence our 
likelihood routine performs an explicit minimization over one nuisance parameter per data 
point. For simplicity, we assume that the error is uncorrelated between different values of 
/, but not between different bins for a given /: this assumption could be relaxed, at the 
expense of increasing the computing time. 

We fixed /max = 2000, since beyond this value both the shot noise term and the 
theoretical error are large, as shown in figure 3. This figure also shows the relative error on 
the observed spectrum in the first and last redshift bins, coming either from observational 
errors (including cosmic variance) or from the theoretical error, and using exactly the same 
conventions as in the previous section: the edges of each of the two error bands correspond 
to a shift between the theory and the observation leading to Ax^ = 1 when either the 
observational or the theoretical error are included in the likelihood. The lowest redshift 
bin incorporates small non-linear scales: this explains why at / = 2000, the theoretical 
error reaches 3.5%. 

Our results are presented in Table 2 for three cases: no theoretical error, uncorrelated 
error only (described in B.4), or additional neutrino-related correlated error (described in 
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Figure 3: Observable cosmic shear power spectrum (top) and its relative error (bottom) for the 
first redshift bin (left) and last redshift bin (right) of a Euclid-X'sk.Q shear survey. The quantity 
displayed above is the lensing auto-correlation spectrum (dimensionless). The upper plots show 
the comparison of a model with massless neutrinos to our fiducial model [My = Zmy = 0.21 eV). 
Solid lines are derived from the non-linear matter power spectrum using the recent update of 
HALOFIT [24], while dotted lines are derived from the linear power spectrum. The lower plots show 
the part of the relative error coming from observational or theoretical errors only (cosmic variance is 
included in the observational error). In these plots, the individual 1-a error on each data point has 
been rescaled by the square root of the number of points, in such a way that the edges of the error 
bands correspond to a shift between theory and observation leading to Ax^ = 1, when only the 
observational or theoretical error is incorporated in the likelihood expression. In these lower plots, 
we also show for comparison the ratio between a massless model and a model with the minimum 
total mass allowed by neutrino experiments. My = 0.05 eV. 



B.5). The impact of the uncorrelated error is again important, but not as pronounced 
as in the galaxy power spectrum case, because on small scales the precision of the shear 
survey is limited by a significant shot noise contribution. The neutrino mass error degrades 
only from a{My) = 0.026 eV to 0.028 eV. For the shear survey we did not perform runs 
with a twice or ten times smaller error: the result for a {My) would simply lie between 
those two numbers. The impact of the neutrino-related error is small but further degrades 
the sensitivity to a{My) = 0.032 eV. While in the absence of theoretical error the galaxy 
survey seems more sensitive to the neutrino mass, the performance of the two methods are 
roughly identical once the same theoretical error ansatz is included. 
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Table 2: Marginalized 1-a error for each raodel pararaeter, in a fit to Planck + Euclid-\ike shear 
survey data. The different lines correspond to the inclusion or not of the global uncorrelated 
theoretical error (un. err.), and of the specific neutrino-related correlated error (co. err.). Our 
preferred prediction is given on the last line, and is very close to that of the second line. 

The triangle plot of figure 4 shows that the parameter degeneracies are very similar 
for the two cases of the galaxy survey and shear survey. Nevertheless, [19] showed that 
combining the two data sets (with a proper cross-correlation matrix) leads to sensitivity 
improvements. It would be interesting to test this conclusion in presence of theoretical 
errors. 

Our result are consistent with those of [11], although a direct comparison is difficult, 
since these authors include several extra parameters (wq^ r, ag) in their forecast. The 
predictions of [19] (case "cs" in their Table 2) he between our results with and without 
theoretical errors. This is consistent since on the one hand, these authors use more op- 
timistic survey characteristics (rf, (t^^^), cTph), and on the other hand, we are including 
much larger values of I (which is legitimate if our theoretical error is realistic). 

4. Conclusions 

We have presented forecasts of cosmological parameters by using, in combination with 
Planck data, two Euclid-like mock future data sets: a galaxy spectroscopic redshift survey 
and a cosmic shear survey. We focused our attention on constraints that can be achieved 
on the total neutrino mass by using the data in the linear and non-linear regimes. 

In order to do this conservatively we adopt the following improvements with respect to 
similar works performed recently in the literature: i) we make use of Markov Chain Monte 
Carlo rather than the Fisher Matrix, which results in more reliable error bars, as well 
as considering degeneracies between parameters. Ultimately, we found that the posterior 
probability is very close to a multivariate Gaussian for the model considered. However, a 
Fisher matrix approach could not have confirmed this, and would not have been explicitly 
independent of the stepsize in the numerical derivatives, ii) we rely on a modification of 
HALOFIT that accounts for massive neutrinos, and predicts the non-linear matter power 
spectrum to small scales, based on the results of N-body and hydro simulations, in) we 
conservatively consider errors both on the non-linear observable power at small scales and 
on the neutrino induced suppression, and explictly show how to implement these errors in 
the likelihood calculation. 

It is instructive to see that with the shape assumed for the uncorrelated theoretical 
error, and a conservative assumption on its amplitude (leading to a 2% error at fcmax = 
0.4/i/Mpc and z = 0.5), the sensitivity to cosmological parameter is still satisfactory. The 
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Figure 4: Marginalized posteriors and two-dimensional probability contours in a fit of Planck + 
Euclid-like shear survey data, with a global uncorrelated error of 5% on non-linear scales (second 
model in Table 2). 



error bar on the total neutrino mass, of the order of 32 meV (cosmic shear) or 25 meV 
(redshift survey), would still allow for a two sigma detection of the total neutrino mass in 
the minimal normal hierarchy scenario. However, with this amplitude and k-dependence 
of the theoretical error, essentially all the information comes from linear scales. The next 
interesting question is to check how much the uncorrelated error should be controlled in 
order to start being sensitive to mildly non-linear scales. Assuming a twice smaller error 
does not change the parameter sensitivity by a significant amount. Extracting significant 
information from non-linear scales requires an error ten times smaller, at the level of 0.2%. 
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Here the error on the neutrino mass decreased from a{M^) — IS meV to 14 meV when 
adding scales with 0.1 < /c < 0.6/i/Mpc to the analysis. This shows that it would be 
extremely useful to be able to predict the observable power spectrum of a given cosmological 
model up to a residual uncorrelated error of the order of 0.1% (resp. 0.2%) at k ^0.1 /i/Mpc 
(resp. k ^0.4/i/Mpc) and z — 0.5. This will be a major challenge for theoretical and 
numerical cosmology in the next decade. 
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A. Galaxy redshift survey implementation 
A.l Observed spectrum 

Let P^^^ be the observed/mock/fiducial power spectrum, and P*^ the spectrum that one 
would expect to see given the theoretical model. Each of these quantities relates to the 
galaxy spectrum Pg and finally to the total non-linear matter spectrum Pnl by taking into 
account redshift distortion effects, spectroscopic redshift errors and light-to-mass bias. A 
good approximation of such a relation is given by (see e.g. [16, 15]): 



DA{z)1^iH{z) th/obs 
Da{zYH{zU ' 



{kref±, kref\\,z) 



(A.l) 



b{zf l + /3{z,k) 




(A.2) 



with the definitions 




)th/obs 
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din a 



{k,z) 



(A.3) 



kre{± = k±H{z)ref/H{z), k^^^ = k\\H {z)re{ / H {z) 
H = kref-r = fcrefll/fcref, 



(A.4) 
(A.5) 




Here b{z) is the bias, assumed to be scale-independent in the range of scales of interest, a 
is the scale factor, H(z) is the Hubble parameter, Da{z) the angular diameter distance, 
and ^{z,k) accounts approximately for redshift space distortions. So we can treat A; as a 



function of the arguments (fcref , /i, z) and write 



3>th/obs / 



^ NL {k(kj;Qf,fi,z),z)e 
A. 2 Likelihood 

For a narrow redshift bin b centered on z, the hkehhood reads 



(A.7) 



A4 exp 
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^ref ref 
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with an effective survey volume given by 



Kff(fcref,M,^) = V;urvey(^) 



ng{z)Pl'^{k^^i,li,z) 
l + ng{z)Pl'^{k^^i,li,z) 



(A.IO) 



survey? ^min 



Later, we wih specify the sensitivity of the survey, parameterized by 
and fcmax- We skip here the derivation of the Fisher matrix, obtained by differentiating the 
above formula twice with respect to the cosmological parameters on which P*^ depends, 
and evaluating this derivative at the maximum likelihood point. We checked that this 
calculation gives exactly the formula commonly used in the literature (see e.g. [16, 15]). For 
the purpose of the discussion in the next section (and also of the numerical implementation), 
we wish to write explicitly the discrete limit of the integrals. We discretize /x in a set of 
equally spaced values /x^, and / = In A: in a set of equally spaced values Ij = In/crefj. The 
step sizes are denoted Afi and Al respectively. We then expand the integral as a sum, and 
for simplicity we omit the factors 1/2 that should weight the boundary terms of each of 
the two integrals. We introduce the short-cut notations: 



Nij = A/xA/ 



^refj^ff (^refj 5 l^ii z) 



(27r2) 



(A.ll) 
(A.12) 



and we get 



Dobs 



(A.13) 



This expression is easy to understand from first principles. Let us consider a single variable 
5 obeying a Gaussian distribution centered on zero and with variance (5^) = P. If we 
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observe N independent realization 6^ of the variable 5, we can build an estimator of the 
variance P of 5, 

(A.14) 

n 

The variance of this estimator can be computed by noticing that each 5^ follows a 
distribution of order one, for which the mean is P and the variance 2P^. So the sum 
J2n^n has a variance 2NP'^. Finally E has a variance {2NP'^)/N'^ = 2P^/N. Moreover, 
E is nearly Gaussian if N is large, as a consequence of the central limit theorem. So the 
probability of the data E given the theory P is a Gaussian of mean P and of variance 
2P^ jN . In other words, 

-21n£(£;|P) = ^^^. (A.15) 

The previous likelihood follows this form for each discrete term. Indeed each term corre- 
sponds to the likelihood of the estimator of the power spectrum in a thin shell in Fourier 
space. The number of independent measurements, i.e. of independent wavenumbers in 
each shell, is given by Niy The role of E and P is played respectively by Pf^^ and P"^. 
Such a likelihood was first derived in pioneering papers like [36, 37]. 

A. 3 Survey specifications 

We computed this likelihood for values of T4urvey(^)7 ^g{^\ ar{z) inspired from currently 
plausible Euclid specifications, which are likely to change over the next years. We divide 
the observations into sixteen redshift bins of width Az = 0.1, ranging from z = 0.5 to 
z = 2.0. For each bin, we assumed: 

• a volume per bin T4urvey(^) — 47r/sky [r(z)]^ (1 + z)~^ ^^^^^ Az, where r{z) is the comov- 
ing distance up to a comoving object with redshift z, with the explicit assumption 
that ao = 1: 

dz' 



H.) = X^. (A.16) 

We assume a sky coverage /sky = 0.375. 

a galaxy number density per comoving volume ng{z), related to the number of galaxies 
per square degree dg{z) through 

4,lr(2)]22gaA.- 

For dg{z), we start from the number presented in Table 2 of [38] for the case of a 
limiting flux of 3 x 10~^^ergs~^cm~^. Following the recommendation of that paper, 
we divide these numbers by 1.37 in order to get conservative predictions. Finally, 
we multiply them by an efficiency factor e = 0.25 (standing for the redshift success 
rate). For instance, for the first redshift bin, this gives dg{z) = 9376/1.37 x 0.25 = 
1710 deg-2. 

a spectroscopic redshift error = ^^§^^z with cr^ = 0.001(1 + z). 
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• a scale-independent linear bias b{z). The choice of b{z) values affects the final result 
less crucially than that of dg{z). We could adopt the predictions of [39] inferred 
from N-body simulations, but for simplicity, our forecast is performed under the 
approximation b{z) = \/l + z. So, we assume in this forecast that the linear bias 
will be accurately measured or predicted for each bin, and that deviations from 
this prediction (coming from the non-linear evolution) will be known up to the level 
described by the theoretical error function. 

• fcmin can be chosen arbitrarily close to zero without changing the results. 

• we tested two values of /cmax^ 0.1 and 0.6/iMpc~^. 

A. 4 Accounting for a global uncorrelated theoretical error 

To present a realistic forecast, one should model all the systematic effects not accounted 
for by the previous likelihood formula, such as: theoretical errors in the calculation of the 
linear and non-linear power spectrum, scale-dependence of the bias on small scales, residual 
shot noise in galaxy counts beyond the contribution already included in the definition of 
V^ff , residual errors in the modeling of redshift space distortion beyond the above scheme. 
On top of these corrections, one may have to take into account the fact that the likelihood 
is not Gaussian on strongly non-linear scales. In this paper, we limit ourselves to mildly 
non-linear scales k < fcmax = 0.6/iMpc~^, and assume that non-Gaussianity effects are 
sub-dominant to the previously mentioned systematics. We also neglect to marginalize 
over residual shot noise in each redshift bin, because Ref. [16, 15] found that this has a 
negligible impact. 

Understanding these various systematics is a major challenge for the future, which 
should be addressed with better simulations and analytical modeling. Here we want to 
keep the analysis simple, and model these systematic errors in a simple way, by adding to 
the spectrum an uncorrelated theoretical error function. By uncorrelated we mean that 
the errors made at different scales are independent from each other, which is the most 
conservative possible assumption. In this case, we can introduce an independent Gaussian- 
distributed nuisance parameter Cij for each data point, and marginalize over it - or rather, 
to a very good approximation, minimize over it: 



where Rij is the theoretical error variance for a bin in (/x, fcref) space centered on (/x^, fcrefj). 
As long as the theoretical error is assumed to be small, it is also a valid approximation to 
neglect the e^j-dependence of the denominator, in order to find a simple analytic solution 
for 6^j, which, injected back in eq. (A.18), gives 




(A.18) 




(A.19) 
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In other words, the theoretical error variance simply adds up to the noise variance. 

Note that we explicitly checked that it is legitimate to neglect the e^j-dependence of the 
likelihood denominator when minimizing over Cij. We also coded the full likelihood with 
explicit minimization over each e^j, and found the same results up to very good accuracy. 

We choose a numerical value of Rij motivated mainly by the current level of precision 
of the HALOFIT algorithm. We assume a relative error on the non-linear power spectrum 
of the form 



a(k^ z) = 



AP^\ik,z) _ ln[l + k/ka{z) 



0.05 



(A.20) 



P^i{k,z) l + \n[l + k/ka{z)] 

where ka{z) is the scale of non-hnearity computed by HALOFIT. This function increases 
from zero to 5% around the scale of non-linearity. Using the function fc(fci-ef , A*, z), this error 
can easily be propagated to the theoretical observable spectrum 



In terms of the discretized observable spectrum, the error reads 

aij = a{k^Qij,iii,z) . 



(A.21) 



(A.22) 



The error variance Rij should be proportional to the power spectrum variance {aijP-'j^)'^. 
We also assume that the error makes a constant contribution to each logarithmic interval 
in the space where observations are performed, i.e. is of the form 



,th\2 



We normalize the error variance Rij in such a way that a one-sigma theoretical error in 
each data point results in increasing the effective by one unit, namely. 
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2B In 
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[aij Pi j 



,th\2 



^refj 



(A.24) 



where B is the number of bins. The role of the normalization factor between squared 
brackets will become clear below. The likelihood becomes (using eq. (A. 19) and going back 
to the continuous limit) 



£ = A/6 exp 



1 d4i 

.1 T 



dk^ei 



ohs 



P 



th\2 



I ^ref^eflP ^min J 



(A.25) 



where we omitted the argument (fcref , ^fe) of the functions P*^, P^^^, Vq{^ and a. If one 
assumes that the observed and theoretical spectra differ by aP^^ for each (A:, /x, z), and that 
in the denominator the theoretical error dominates over the observational one (V^ff = oc), 
then 



£ = A4 exp 



_i dfi 



dkj-ei 



1 



Kmin 



(UbMb) exp 



(A.26) 
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which corresponds to a shift by AXgf^ = 1 with respect to the maximum hkehhood L — 

If we had assumed the error to be fully correlated, instead of increasing the denominator 
of the likelihood, we would have replaced P*^ by P*^(l + ea), multiplied the likelihood by 
072^exp[-eV2], and marginalized/minimized over e. Then, the assumption P^^^ — 
P*^(l + a) would correspond to an optimal choice e = 1 in the large V^ff limit, and would 
also lead to a shift in Ax^f^ by one unit with respect to the assumption P^^^ = P*^. In our 
case, we obtain the same shifting while assuming statistically independent errors for each 
data point. 

Finally, the likelihood can be simplified to 



C II5A4 exp 



ref 



fcref (27r)2 



ref "^survey 



Aref 



'obs 



H pth 



JL pth 
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D\ rig 
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Pin 

'^mir 

(A.27) 



where we omitted the argument z^, in the functions V^urvey, Da^ H and Ug. This is exactly 
the relation implemented in our code. 

A. 5 Accounting for an extra neutrino-related error 

The impact of massive neutrinos on non-linear corrections to the power spectrum has been 
investigated in [24]. By comparing with N-body simulations including neutrino particles, 
the authors of [24] re-calibrated HALOFIT, with a new neutrino mass dependent correction. 
This fitting procedure is of course not perfect and adds a systematic error growing with 
the neutrino mass. It was found that the leading error can be described with a correction 



Hi + k/kAz)] 



(A.28) 



with fi, = ijjv/ojm-, and e^, is an unknown correction of unit variance, that we will treat as 
a Gaussian nuisance parameter. Hence our final definition of the likelihood accounting for 
both types of error reads 



C — M exp 



(P°'^^-[P*ni + e,(7,)])^ 



)th 



exp 



A^ref^eff ^min 



(A.29) 



where we omitted the argument (fcref , A^, z) of the functions P^^^, P*^, and V^fT- Note 

that the correction proportional to should not be added to P°^^ since we are assuming 
for simplicity that the fiducial value of in the mock data is zero. 
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B. Cosmic shear survey implementation 



B.l Observed spectrum 

As in e.g. [16, 15], we define the likelihood of the shear auto or cross-correlation power 
spectrum in bins i and j: 

Here, Wi{z) is the window function of the z'th bin. It can be evaluated as a function of 
the radial distribution of galaxies in each redshift bin, Di(z), obtained by convolving the 
full radial distribution D{z) with the photometric redshift uncertainty function P(^,^ph), 
multiplied the top-hat window function of each bin: 

Wi{z) = 5si„(l + z)Fi{z) (B.2) 

F.(.)= r!MM£MzlM,,, ,B.3) 

Jo r[Zs) 

Jo Di{z')dz' 

D,{z) = / ^ P(2:,2:ph)i?(zph)rf2;ph . (B.5) 

The radial distribution D(z) can be arbitrarily normalized, since ni{z) is anyway normalized 
to one. We will assume that the photometric redshift uncertainty function is normalized 
/o^ ^(^5 ^ph)^^ph — I5 but a different normalization would not impact the final result 
for the same reason as for D(z). The noise spectrum contaminating the measurement of 
Ci"^ is given by the diagonal matrix in ij space: 

Nr^S,,{j^,^,)n-\ (B.6) 

where {"^rms)^^^ ^^^^ mean square intrinsic shear (like in the forecasts of the Euclid 

Red Book [1], we assume that this quantity is equal to 0.30), and Ui is the number of 
galaxies per steradian in the i'th bin, given by 

Hi = 3600d(180/^)^ni , (B.7) 

where d is the full number of galaxies per square arcminute in all bins, and hi is the fraction 
of galaxies in the z'th bin, given by: 

J min ^(^) 

n^ = , • (B.8) 

Jo ^(^) 

We used the survey specifications for D(z)^ d and fsky detailed in Appendix B.3. 

Using dz/dr — i7, we can write the same integrals in a different way (used in other 
papers and in our code): 

C]^ = ^^^iHl drr-^g.{r)g,{r)p{k = ^,^(r)) (B.9) 
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with 

g,{T) = 2r(l + z{r)) T ^^^ ^^(^^)(^^ " ^) (b.10) 
rji{r) = H{r)n,{z{r)) (B.ll) 

and ni{z) is the same as before. 

B.2 Likelihood 

Let's assume some theoretical spectra Cf^^^ (here, the spectra of each model that we want 
to fit to the data, exploring the space of free cosmological parameters), and some observed 
spectra C^^^^K The matrix C° ^ of element C^^^^^ is called the data covariance matrix. It 
can be inferred from the observed multipoles a^}^'*'^ which are Gaussian distributed with a 
variance independent of m in an ideal full-sky experiment, so that 

Cf^'^ = {21 + Yl • (B.12) 

m=—l 

For a parameter forecast, instead of the covariance matrix of mock data, we can use some 
fiducial spectra corrected by the noise spectra of the experiment at hand: 

^obsij _ ^fiducialij ^ ^g^^g^ 

This data covariance matrix should be compared with the theoretical covariance matrix 
defined as 

^thij ^^thij^^ZJ ^^^^^ 

We define the determinant of these N x N symmetric matrices: 

df" = det (C^^'^^ (B.15) 

df' = det(c;^^'^j) . (B.16) 

The determinants are homogeneous polynomials of order N in the spectra, e.g. for N = 2: 

df = cf^^^cf'^'^ _ ((5*^12)2 _ (B.17) 

The quantity df^^^ can be built starting from df^^ and replacing one after each other the 
theoretical spectra by the corresponding e.g. for N = 2: 

^mix ^obsll^th22 _|_ ^thll^obs22 2^thl2^obsl2 /g jg") 

I I I I I II \ J 

So, df^^^ is always linear in the C^^^^^^s. By construction, when C^^^^ = C^^^^^ one has 
^mix _ ^^th _ TYd^^^. Since in an ideal full-sky experiment, the different multipoles are 
uncorrelated in (/,m) space, the likelihood of the observed spectra given the theoretical 
spectra is as simple as: 



1 obst//=ith\-l obs 



(B.19) 
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where a^^^ = {^^m^} N-dimensional vector of observed multipoles in each bin, Cj^^ is 

the theoretical covariance matrix of element and A/" is a normalisation factor. After 

some simple algebra^, the likelihood simplifies to 



The effective chi square 



\ (41^)1/2 



exp 



{21 + 1) 
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reaches its minimum for C°^* — C*'*, corresponding to 

Xlr^" = -2 In £„ax = -2 In + 5^(2^ + 1) (iV + In df') 
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The relative to the best-fit model is then equal to 
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(B.20) 



(B.21) 



(B.22) 



(B.23) 



equivalent to 
(B.24) 



Finally, a first-order approximation to account for the limited sky coverage of a given 
experiment, consists of increasing the cosmic variance by a factor /s^y^^ 
postulating: 

/ jmix Jth 

AXe\ = B2^ + l)/sk.(^ + lnA_-iV 

This is precisely the expression used in the code. 
B.3 Survey specifications 

A given survey is specified by D(z), and finally by the covered faction of the sky 

fsky'i it can then be decomposed in redshift bins according to some strategy defined by the 
user. For a Euclid-like experiment we use the same characteristics as in the Euclid Red 
Book [1]: 

exp[-(z/zo)^-^] for z < z^^^ = 3.5 
with mean redshift Zn 



Diz) 
'P{z, ^ph) 
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fsky 



1 



: exp 



-mean - IAUzq = 0.9 
21 



(B.25) 
(B.26) 



with cTph = 0.05(1 + z) 
30 arcmn 
0.375 . 



(B.27) 
(B.28) 



We assume five bins, with the first bin starting at z^^^^ = 0, the last one ending at z 
3.5, and bin edges z^^^^ = z^^ 
galaxies, i.e. fii = 



max 

N — 

chosen such that each bin contains the same number of 



^in particular, using A = adj(A)/ det{A) where adj(A) is the adjugate matrix of A, i.e. the transpose 
of the matrix of cofactors of A. 
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B.4 Accounting for a global uncorrelated theoretical error 

Like for the power spectrum likelihood, taking into account an uncorrelated error on each 
data point is equivalent to minimizing over a number L = (/max — /min + 1) of nuisance 
parameters e/: 



-cx)<e<+cx) 

/ — ^min 



(B.29) 



Here, $j^{ei) stands for the determinant of the theory covariance matrix shifted by the 



theoretical error covariance matrix : 

Jf(Q) = det(C;^'J + qi?n. (B.30) 

Similarily, (ip^^(e/) stands for the sum of N terms, each one being the determinant of a 
matrix built from ^■^ + ^iR\\ where one column has been replaced by the same column 
in the observed covariance matrix. Hence the quantity d™^ defined just above eq. (B.18) 
is identical to d^'^'iO). 

Note that for simplicity, we consider here uncorrelated errors for each Z, but not for 
each bin. This approach could easily be generalized to independent bin errors, at the 
expense of introducing more nuisance parameters. 

In the case of the power spectrum likelihood, we could find an analytical approximation 
of the nuisance parameter value minimizing the effective . In the present case, we checked 
that simple approximate solutions are not accurate enough. We perform a numerical 
minimization over each ei within the likelihood routine, using Newton's method. 

We define our theoretical error covariance matrix i?^-^ in a similar way as for the power 
spectrum likelihood. We start from the power spectrum relative error function a{k^z) 
defined in eq.(A.20). The power spectrum error can be propagated to a covariance matrix 
error E^j^ \ 

E'/ = ^nlH^ J^^ drr-'gi{r)gj{r)al^k = ^-,z{r)^ P'"^ (fc = ^-,z{r)^ . (B.31) 

The theoretical error matrix i?^-^ should be proportional to E^^ . We normalize it to 

H^^^ ^L^I^E]^ (B.32) 

in such a way that enforcing a one-sigma theoretical error for each I results in an increase 
of the by one (as would be the case for a fully correlated theoretical error with the 
same amplitude). Then, if one assumes that for each / the observed spectra are equal 
to the theoretical ones shifted by a one-sigma theoretical error [Cf^^ — Cf^ + ^r)? 
minimization gives (up to a very good approximation) ei — L~^/^, and 

^Xeff = E [(2^ + l)/«ky (iV + - iV) + = 1. (B.33) 
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B.5 Accounting for an extra neutrino-related error 

Finally, we account for the correlated error modelling neutrino-related uncertainties by mul- 
tiplying the theoretical power spectrum P*^(fc, z) by a factor (l + ejyajy(A;, z)), as in equation 
(A. 28), as well as adding to ^xlfi- The nuisance parameter e^y is then marginalized over. 
Note that the factor (1 + ejya^{k^z)) should not multiply the observed/fiducial spectrum, 
as long as we assume a fiducial value of e^y equal to zero. 
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